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^~-) ■ Abstract. We study the Einstein- Vlasov system coupled to a nonlin- 

, ear scalar field with a nonnegative potential in locally spatially homoge- 

neous spacetime, as an expanding cosmological model. It is shown that 
solutions of this system exist globally in time. When the potential of 
the scalar field is of an exponential form, the cosmological model corre- 
sponds to accelerated expansion. The Einstein- Vlasov system coupled to 
a nonlinear scalar field whose potential is of an exponential form shows 
the causal geodesic completeness of the spacetime towards the future. 
, The asymptotic behaviour of solutions of this system in the future time 

is analysed in various aspects, which shows power-law expansion. 

o ' 
^t- 
o 

O ■ 1. Introduction 

O 

^ Particle systems are modeled statistically by distribution functions, which 

^jTj| at any time represent the probability to find a particle in a given position, 

with a given momentum. The distribution functions contain a wealth of 

information and macroscopic quantities are calculated from these functions. 

The models being considered here are those in which collisions between par- 
es ' 

tides are sufficiently rare to be neglected. The collection of these collisionless 
particles is described by Vlasov equations. For this reason, matter consid- 
ered in these physical models is said to be collisionless matter or Vlasov 
matter. 

The time evolutions of particle systems are determined by the interactions 
between the particles which rely on the physical situation. Each particle is 
driven by self-induced fields which are generated by all particles together. 
Naturally combinations of interaction processes are also considered but in 
many situations, one of them is strongly dominating and the weaker pro- 
cesses are neglected. In gravitational physics, these fields are described by 
the Einstein equations. The physical models concerned in this paper is de- 
scribed by the Vlasov equation which is coupled to the Einstein equations 
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by means of the energy-momentum tensor. One application of the Vlasov 
equation coupled to this self-gravitating system is cosmology. The particles 
are in this case galaxies or even clusters of galaxies. 

The simplest cosmological models are those which are spatially homoge- 
neous. Spatially homogeneous spacetimes can be classified into two types; 
Bianchi models and the Kantowski-Sachs models. The models with a three- 
dimensional group of isometries G3 acting simply transitively on spacelike 
hypersurfaces are Bianchi models. There are nine types I-IX, depending 
on the classification of the structure constants of the Lie algebra of G3. 
Those admitting a group of isometries G4 which acts on spacelike hypersur- 
faces but no subgroup G3 which acts transitively on the hypersurface are 
Kantowski-Sachs models. In fact, G3 subgroup acts multiply transitively on 
two-dimensional spherically symmetric surfaces. 

If we take as a cosmological spacetime one which admits a compact 
Cauchy hypersurface, the Bianchi types which can occur for a spatially ho- 
mogeneous cosmological model are only type I and IX and also Kantowski- 
Sachs models. Because of the existence of locally spatially homogeneous cos- 
mologies, we take a larger class of spacetimes possessing a compact Cauchy 
hypersurface so that this allows a much bigger class of Bianchi types to be 
included. Since the Cauchy problem for the Einstein- Vlasov system is well- 
posed, it is enough to define the class of initial data. Here is the definition. 

Definition 1. Let g^, kij and JT be initial data for a Riemannian metric, a 
second fundamental form, and a matter, respectively, on a three-dimensional 
manifold M. Then this initial data set (9ij, kij, J 7 ) for the Einstein- Vlasov 
system is called locally spatially homogeneous if the naturally associated 
data set on the universal covering M is homogeneous, i.e., invariant under 
a transitive group action. 

So the spacetimes considered here will be Cauchy developments of locally 
homogeneous initial data sets on some manifolds. Note that a complete 
Riemannian manifold is locally homogeneous if and only if the universal 
cover is homogeneous. For Bianchi models the universal covering space can 
be identified with a Lie group G. So the natural choice for G in this case is 
a simply connected three-dimensional Lie group. (For a detailed discussion 
on this subject we refer to Ej)- 
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In this paper, we discuss the dynamics of expanding cosmological mod- 
els, particularly accelerated expansion. There are two subjects concerning 
this rapid expansion. One is the very early universe close to the big bang 
(inflation) and the other is the present era (quintessence) supported by the 
observations of supernovae of type la. 

One simple way to obtain accelerated expansion is to introduce a positive 
cosmological constant, which leads to exponential expansion. In homoge- 
neous spacetimes it has been studied by Wald in jllj with general matter 
which satisfies the dominant and strong energy conditions. When the matter 
is described by the Vlasov equation, the detailed asymptotics of solutions 
have been analysed in [I]. In the inhomogeneous case Vlasov matter model 
has been studied in under some symmetric conditions. In [§] by Ren- 
dall, vacuum and perfect fluid cases are handled. 

Another choice for accelerated expanding cosmological models, which is 
more sophisticated, is a nonlinear scalar field. It has been analysed by 
Rendall in ^01 that when the potential of the scalar field has a positive 
lower bound with general matter satisfying the dominant and strong energy 
conditions then the homogeneous models expand exponentially. In the case 
of an exponential potential, the models shows power-law expansion which 
has been studied in [21 0] by Kitada and Maeda. 

Bianchi type IX and Kantowski-Sachs models have complicated features 
when a positive cosmological constant or a nonlinear scalar field is present. 
It has been seen that there are chaotic behaviours between expanding and 
recollapsing phases in these models. In the discussion of expanding cosmol- 
ogy, the models being concerned in this paper are all Bianchi types except 
IX, with a nonlinear scalar field and vanishing cosmological constant. 

The fields described by the Einstein equations are coupled to the Vlasov 
equation and a nonlinear scalar field by the energy-momentum tensor which 
is of form 

T ^ = T^ lasov) + V a cf>V^ - (iv^V 7 + V{<j>))g a p (1.1) 

where <f> is a scalar field which represents dark energy, V is a potential 
and j"jy iasov ^ [ s th e energy-momentum tensor of the collisionless matter de- 
scribed by the Vlasov equation. r p < y^ asov ^> satisfies the dominant and strong 
energy conditions given respectively by 
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(1) T^a as v a w^ > where v a and are any two future pointing timelike 
vectors, 

(2) (T^ lasov) - y a p(T$ lasov) gn)v a vP > for any timelike vector v a 

As a consequence of the Bianchi identity in the scalar field <j> satisfies 
the equation 

V a V a <P = V'(</>) 

To make the notation not too heavy the superscript (Vlasov) will be omitted 
for the rest of the paper. 

The content of the rest of this paper is the following. Section|2]presents the 
detailed formulation of the system being considered. In Section we prove 
the global existence of solutions for the Einstein- Vlasov system coupled to 
a nonlinear scalar field with a general potential. Also the causal geodesic 
completeness of the spacetime towards the future will be present, in the 
case of an exponential potential. In Section 21 we study the asymptotic 
behaviour of solutions at late times in various aspects, when the potential of 
the scalar field is of an exponential form. We observe the future asymptotic 
behaviours of the mean curvature, the metric, the momenta of particles 
along the characteristic curves as well as the generalized Kasner exponents 
and the deceleration parameter. Also we analyse the energy-momentum 
tensor in an orthonormal frame on the hyper surf aces. As we will see later 
on, the cosmological model being considered in this paper exhibits power-law 
expansion. 

2. Einstein- Vlasov system with a scalar field 

Here is the formulation of the Einstein- Vlasov system coupled to a non- 
linear scalar field with a potential. Let G be a simply connected three- 
dimensional Lie group and {e^}, a left invariant frame and {e*}, the dual 
coframe. Consider the spacetime as a manifold G x I, where I is an open 
interval and the spacetime metric of our model has the form 

ds 2 = -dt 2 + gijity ® e j (2.1) 

The initial value problem for the Einstein- Vlasov system is investigated in 
the case of this special form of the metric and the distribution function / 
depends only on t and v l , where v % are spatial components of the momentum 
in the frame e,. Initial data will be given on the hypersurface G x {to}- 
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Now the constraints are 



R - (kijk ij ) + {k ijg ij ) 2 = 167rT o + 8vn/> 2 + 16ttV(<£) 
V*fcy = -8nT j 



(2.2) 
(2.3) 



The evolution equations are 



d 




(2.4) 




Rij + (ki m g lm )k tj - 2k u k] - wry 

- 47rT oogij + 4ir(T lm g lm ) gi] - SirV^g,, 



(2.5) 




(k lm g lm )*P ~ V\<t>) 



(2.6) 



(2.7) 



These equations are written using frame components. Here kij is the second 
fundamental form, R is the Ricci scalar curvature and Rij is the Ricci tensor 
of the three-dimensional metric. And (ft is a, scalar field, depending only on 
t, with a nonnegative potential V((ft). tft is a function introduced by the 
relation (|2.6|) . The nonnegative assumption on the potential is very natural. 
It implies that the dominant energy condition is satisfied and then the weak 
energy condition follows. 

Here are the components of the energy-momentum tensor of the Vlasov 
matter ; 




1 +g rs v r v s ) 






(2.10) 



(2.9) 



Here v := (f 1 ,^ 2 ,^ 3 ) and dv := dv 1 dv 2 dv s . 
The Vlasov equation is 



d t f+{2kjvi - (l+g rs v r v s r 1/2 f mn v m v n }d vl f = 



(2.11) 



Here the Ricci rotation coefficients r f mn are defined as 
7mn = 7;9 tk (—C l n k9ml + Cfcm^ni + C l mn , 
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where Cj fc are the structure constants of the Lie algebra of G. To have a 
complete set of equations it is necessary to compute Rij in terms of gij. In 
this paper, it is enough to know that is of the form (det g)~ n (polynomial 
in gij and Cj fc ). To control (detg) -1 we use 

^log(det 5 ) = -2(%^) (2.12) 

Note that in discussing expanding cosmological models, the sign conven- 
tion for {kijg 13 ) in the paper is negative. Also it is true that if models are 
initially expanding, i.e., (kijg lJ )(to) < then (kijg l3 )(t) < for all time 
t > t (For details see [TU]l 

The evolution equations are in general partial differential equations, i.e. 
d/dt is dt- However due to the locally spatially homogeneous spacetime, the 
partial differential equations are reduced to ordinary differential equations. 

For the rest of the paper, C denotes a positive constant which changes 
from line to line and may depend only on the initial data. Also C\ (I = 
0, 1,2,.. .) are positive constants. 

3. Global existence of solutions and geodesic completeness 

In this section, we will show global existence solutions of the Einstein- 
Vlasov System coupled to a nonlinear scalar field with a potential. As a first 
step, conditions will be established under which solutions of this system exist 
globally in time, with the technique appeared in [Jj by which the existence 
of solutions for the Einstein- Vlasov system in the absence of a scalar field 
has been proved. And then eventually it will be proved that these conditions 
are fulfilled in the system being considered. Also we will observe the casual 
geodesic completeness of the spacetime towards the future direction, when 
the potential of the scalar field is of an exponential form. 

Proposition 1. Let g^ (to), kij(to), <f>(to), VK^o) and f(to,v) be an initial 
data set for the evolution equations \2.J$ - \2. 7p and the Vlasov equation 
\2.11\) which has Bianchi symmetry and satisfies the constraints \2.2\) and 
\2. c J\j) . Also let f(to,v) be a nonnegative C 1 function with compact support. 
And assume that the potential of the scalar field V(4>) is a nonnegative C 2 
function. Then there exists a unique C 1 solution (gij, kij, (j), ip, f) of the 
Einstein- Vlasov system, on an interval [to,T), for some time T. If \g\, 
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(detg) , \k\, \4>\, \tp\, \\f\\ and the diameter of supp f are bounded on 
[Iq, T), then T = oo. 

Proof : The characteristics of (|2.1ip are the solutions V l (s,t,v) of the 
equation 

A\ri 

— = 2k)Vi - (1 + g rs V r V s r 1/2 lLv m V n (3.1) 

with V%t,v) = v\ Let f<f>\t,v) = f(to,v), g§\t) = 9ij (t ), k§\t) = 

kijfa), = 0(t o ) and V (0) (*) = V>(*o). If / (n) , g§\ *§\ <A (n) and 

are given for some n, determine by solving the characteristic 

equation with k$ and g§\ Let f {n+1 \t,v) = /(*o,V( n+1 )(* , *,«))■ 

Define an energy-momentum tensor T^~^ with y( n + 1 ) and gffi in (|2.8|) - 
(t2~TTH) . Determine g% +1) , k^ +l) , and <i/? (n+1) by solving (J23J) - (t2~71) 

with , <7y , fey 5 0^ an d m the right hand side of equations 

and with c/J" +1) , fe-" +1) , c/>( n+1 ) and ^ (n+1) in the left hand side. Now let 
[io,T( n+1 )) be the maximal interval on which ffjj is positive definite. By 
induction, one can see that /(»), ^ n) , fcg°, and V W 

are C 1 on their 

domains of definition. 

Let \g\ be the maximum modulus of any component g^ and |fc| for fey. 
Suppose that for all n < N — 1 the following bounds hold: 

\g {n) -<7 (0) l <Ai, (det^)- 1 < A a , |fe^-fe^|<A 3 (3.2) 
_ 0(0) | < Aai \^{n) _ ^(o)| < Ag (3 _ 3) 

Also suppose that |u| < A$ whenever f^ n \t,v) ^ 0. Here A{ (i = 1, ... ,6) 
are positive constants which are for the moment arbitrary. The character- 
istic system ()3.1() implies a bound for the form \v\ < Cq + whenever 
f( N \t,v) ^ 0, where i?6 depends only on ij's. As a consequence (|2,8jl 
- (|2.10|) imply a bound for T^jj depending only on A^s. The evolution 
equations 1)2 .4|) - (|2.7j) imply bounds of the form 

\gW-gW\< Bl t, |feW-fe(°)| <B 3 t 

|0W-^°)|<fl 4 t, |V W -V (0) I <^ 5 t 

where Sj's depend only on Aj's. If Ai's are fixed then the inequalities in 
(|3.2() imply an inequality of the form (det g^-*) -1 < -B2 whenever t <T and 
T is some positive time depending only on Ai's. Now fix A^s in such a way 
that A2 > (det i/ - 1 ) -1 and ^6 > Cq. Next reduce the size of T if necessary 



HAYOUNG LEE 



so that BiT < Ai (i=l, 3, 4, 5), B 2 < A 2 and C + B 6 T < A 6 . Then all 
iterates exist on the interval [t ,T) and g( n \ k^ n \ 4> {n) and V (n) are bounded 
on that interval independently of n. 

Now we need to show that these iterations converge. Consider the differ- 
ence of successive iterates for n > 2, 

\(g (n+1) -g {n) )(t)\ + \(k^-k^)(t)\ 

+ |(^ n+1 ) - <t>W)(t)\ + |(0 (n+1) - ^ (n) )(t)l 

)(g {n) -9 (n ~ 1) )(s)\ + \(k^-k^)(s)\ + 



< c 



to 



+ \& {n) - ^ n - 1] )(s)\ + ll(/ (n+1) " f in) )(s)\\ 
For the difference of the characteristics note that 



ds 



(3.4) 



_ jLy(n) < criy(n+l) _ y(n)| + \ g (n) _ g (n-l) , + i fc (n) _ fc (n-l) H 

ds <is 

(3.5) 
(3.6) 



Define 

a( n >(*) := \(g {n+1) - g {n) )(t)\ + |(fc (n+1) - k^)(t)\ 
+ |(<^ +1 ) - (n) )(t)| + |(^ (n+1) - ^ (n) )WI 

+ sup{|y( n+1 ) - V^\(s,t,v) : s G [t ,t], u G supp/(" +1 )(t) U supp/ <%)} 
Then we get 

||/ (n+1) (*)-/ (n) (t)||oo<||/ (0) || O iaW(t) 
Therefore - (E3) imply that 

a (n) (t) <C f [ 

Jto 

Applying Gronwall's inequality to this gives 



(3.7) 



|aW(i) + oM( s )] (i s 



a (n) (*) < C 



a 



(n-l), 



s) ds 



Therefore a( n )(t) < C n ~ 2 \\a^ \\t n ~ 2 /(n - 2)1 and so {g {n) }, {k^}, {4> {n) }, 
{^ (n) } and {VW} are Cauchy sequences on the time interval [to , T) . Denote 
the limits of these sequences by g, k, <f>, ip and V^, respectively. Also by (|2.4|) 
- (|277|) . dg^/dt, dk^/dt, d^/dt, dip^/dt and dV^/dt are uniformly 
convergent. Thus (g, A;, 0, ■0, /) is a C 1 solution of the system on the interval 
[to,T). 
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Now let us check whether a solution exists uniquely or not. If two solutions 
with the same initial data are given, define a quantity a(t) in terms of their 
difference in the same way that aW(i) was defined in terms of the differences 
of two iterates. Applying the same argument as above leads to an estimate 
of the form 

a(t) <C a(s)ds 
Jt 

By Gronwall's inequality we can see that a(t) is zero and hence that the 
two solutions agree. Therefore the solution which has been constructed is 
uniquely determined by the initial data. 

Define 

A := R - {hjk ij ) + (k ijg ij ) 2 - levrToo - 8vn/> 2 - 167rV(^) 
Ai := V%- + 87rT 0j 
Then after a lengthy calculation we obtain 

j t A = 2(k l3 g^)A-2 1 \^A l 

j t Ai = (k lm g lm )Ai + 2k\At 

That is, ()2.4|) - (|2.11|) imply a homogeneous first order ordinary differential 
system for constraints (|2.2|) and (|2.3|) . Therefore we can conclude that if the 
initial data satisfy the constraints then so does the solution of the evolution 
equations (|2.4j) - (|2,7|) with energy-momentum tensors (|2.8|) - (j2.1U|) and 
the Vlasov equation (|2.11|) . 

In the above argument so far, we see that the size of T is only restricted 
by the quantities ; \g(% (det^ ))" 1 , |Jfe(°)|, \<j>i% |^ {0) |, ||/ (0) || and the 
diameter of supp f(°\ Thus if the quantities \g\, (det^) -1 , |fc|, \<j)\, \\f\\ 
and the diameter of supp / are bounded on the same time interval [t ,T), 
then a solution exists on [t, t + e) for any t S [to, T) and some e independent 
of t. It can be concluded that the original solution can be extended to the 
larger interval [to,T + e). Therefore this completes the proof. □ 

Let us state some properties of linear algebra which can be found in j^|7j. 
We shall make use of these properties later on. Let ibeanxn matrix. Let 
Ai and A2 be n x n symmetric matrices with A\ positive definite. Define a 
norm of a matrix by 

\\A\\ := sup{Px||/||x|| : x / 0, x G K™} 
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Also define a relative norm by 

||A 2 |U a := sup{||A 2 s||/||Aia;|| :^0,xel"} 
Then from these definitions, one can see that 

||A 2 || < ||A 2 |Ui||Ai|| (3.8) 

and also 

\\A 2 \\ Al < (tr(A^A 2 A^A 2 )) l/2 (3.9) 
Proposition 2. If (kij g^) is bounded on [to,T), then T = 00. 

Proof : Let be the trace free part of the second fundamental form 
kij. Then we have kij = ^(kijg^)gij + &ij. By this fact, we rewrite the 
constraint (12.21) as 



\(hj9 ij f = -\R + + SttToo + 4ttV 2 + 8irV(<f>) (3.10) 

It has been proved by Wald in 111] that in all Bianchi models except type IX, 
the Ricci scalar curvature is zero or negative. Also due to the nonnegative 
potential, we get 

So if (hjg 11 ) is bounded on [to, T), then ip is bounded on [to, T) and so is <fi. 
From evolution equations (|2.4j) and (|2,5|) . we have 

j t {h 3 g ij ) = R + {ki 3 g ij ? + IrSY},;/') - 12nT 00 - 24^(0) (3.11) 

Using the constraint (|2.2[) we get 

d 
dt 

Thus we obtain 



(k ij9 ^) = (h 3 k ij ) + MTij9 ij ) + 47tT o + 8tt^ 2 - 8nV(<f>) 



Then 



(%^)+8tt / V((f>)(s)ds>k ij (t )g l ^t )+ / (kijV3)(s)ds 

J to J tg 



Note that 



/ V(cj))(s)ds <\\V'\\ t [ \<j>(s)\ds + C(t + l) 

Jto J t 
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where [|V||t := sup{| V(</>)(s)| : for all s E [io>*]}- Since 4> ls bounded on 
[to, T), then V(<p)(s) ds is bounded for all t in [to, T). Therefore with the 
boundedness of (kijg 1 ^), we conclude that \kijk l i)(s) ds < 



oo. 



Let \\g\\ and ||fc|| be the norms of the matrices with entries g^ and kij, 
respectively. Let \\k\\ g be the relative norm of the matrix with entries kij 
with respect to the matrix with entries g^. Then using (J3.9|) we have 



< \\g(t )\\+2 fields 

J t 

< \\g(to)\\+2 [ \\k(8)\\ g \\g(8)\\da 

J to 



< ||^o)||+2 / (%^') 1 / 2 ( S )|| 5 (s)||d S 
By Gronwall's inequality, we get 

|| fl (i)|| < ||<?(*o)||exp [2 f(k ij 0) l l 2 (s)ds 



L J to 

Since J* (kijk^) l ^ 2 (s) ds is bounded on [to, T), also \g\ is bounded on [to,T). 
Using (|2.12|) . we see that (detg) -1 is bounded on the same interval. It is 
known that if (det g) and its inverse are bounded then the scalar curvature 
R is bounded from above. Note that in (|2.2I) we have 

R+ihrfif > (hjkV) 

Thus kijk^ is bounded on [to,T). By the inequality 

||fc|| < (%^') 1/2 ||<?ll 

also |fc| is bounded. The boundedness of \g\ and (detg) -1 implies that g 
is uniformly positive definite on the interval. Hence the solutions of the 
characteristic equation are also bounded. Therefore by Proposition the 
proof completes. □ 

So far it has been proved that solutions of the system exist as long as 
some quantities are bounded in a finite time interval [to,T) for arbitrary T. 
These conditions are satisfied, as we will see in the following theorem. 

Theorem 1. Let gij (to), kij (to), 4>(to), ip(to) and f(to,v) be an initial data 
set for the evolution equations JlO ) - and the Vlasov equation h2.ll]) 

which has Bianchi symmetry and satisfies the constraints h2.2]) and h2.3\) . 
Also let f(to,v) be a nonnegative C 1 function with compact support. And 
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assume that the potential of the scalar field V{4>) is a nonnegative C 2 func- 
tion. Then there exists a unique C 1 solution (gij , kij ,<p,tjj,f) of the Einstein- 
Vlasov system for all time. 

PROOF : Consider ^TW\ with (1X1171) 

J t ( k v9 ij ) = ~R + IfajV*') + MTij9 ij ) + 127tToo + 12vr^ 2 

Then 

J t (hj9 ij ) > (3.12) 

Since the cosmological models we are considering here is expanding, i.e., 
{kijg 13 ) < 0, with (|3.12|) we conclude that (k^g 13 ) is bounded for t > to and 
the proof completes from Proposition |2 □ 

3.1. Geodesic completeness with an exponential potential. The next 
result asserts the geodesic completeness of locally spatially homogeneous 
spacetimes for the Einstein- Vlasov system coupled to a nonlinear scalar field 
whose potential is an exponential form. 

Theorem 2. Suppose the hypotheses of Theorem^ hold. And assume that 
the potential of the scalar field V((f>) is of form 

V(<f>) = V e- XK * 

where Vq is a positive constant, A G (0, y/2) and k 2 = 8tt. Then the spacetime 
is future complete. 

The proof of this theorem can be founded in Subsection 14.71 

4. ASYMPTOTICS OF SOLUTIONS WITH AN EXPONENTIAL POTENTIAL 

We study the asymptotic behaviour of solutions in the future time with 
a particular form of the potential V ((/>)■ Namely the potential is given, as 
in the previous section, by V{<p) = Voe~ Aft< ^ where Vq is a positive constant, 
< A < v2 and k 2 = Sir. In has been shown in that in order for power- 
law inflation to occur A must be smaller than v2- Note that the case A = 
corresponds to the model with a positive cosmological constant instead of 
the scalar field which has been well understood in [I]. Briefly, this model 
exhibits exponential expansion. For detailed information, we refer to [3]. 
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We introduce a new time coordinate r denned by 

dr = e~ XK(t>l2 dt (4.1) 

And let kj ■= %e AK< ^ 2 , R := Re x ^ ', f a p := T Q/3 e A ^ and ^ := ^e AK< ^ 2 . 
Then the Hamiltonian constraint (12.21) become 



R - {hjk ij ) + (k ijg ij f = levrToo + 8vrV> 2 + 16irV (4.2) 
The evolution equations are 

^-gij = -2kij (4.3) 



dr 
d 

7h 



kj = Rij + {hm9 lm )~hj - 2kak\ ~ 87rTy (4.4) 
- 47r?oo#ij + 4:TT(fi m g lm )gij - SnVogij + ^-kjip 

±.<j> = & (4.5) 

= ( + AkV + (4.6) 
Also the Vlasov equation becomes 

d T f + {2ky - e A ^ 2 (l + g rs v r v s r 1/2 YmnV m v n }d v ,f = (4.7) 

Now we define two functions : 

_ _ . . 2 - 

S(^,kijg ij ) ■= (hg lj ? - 12ir(^ 2 + 2F ) 

Note that the function S will play the same roles as (ktjg ZJ ± 3A) in the 
papers [H ITT]. 

The basic idea of the following proposition is from [3] . Here the computa- 
tion is carried out carefully so that the error terms are explicitly determined 
for the future reference. 

Proposition 3. Let <Jij be the trace free part of the second fundamental 
form kij such that 

kj = ^(k m g lm )gij + Oij (4.8) 



14 HAYOUNG LEE 



where a^j := aije XK ^^ 2 Then we have 



S = <D{eT e T ) (4.9) 

= 0(e- £ * T ) (4.10) 

R = 0{e- € * T ) (4.11) 

Too = 0{e~^ T ) (4.12) 

Proof : Note that using (j4.8|) we rewrite the constraint (|4.2|) as 

{hj9 13 ? = ~\ R + \{*ij* iS ) + 247r:f oo + 12vrV; 2 + 24ttV (4.13) 
So S becomes 

S$, kij9 ij ) = -\ R + + 247rToo (4.14) 

Recall that the Ricci scalar curvature is zero or negative in all Bianchi models 
except type IX. So the models are allowed to evolve in the region of S > 0. 
Note that from (|4.3j) and (|4.4j) we have 

^{hi9 ij ) =R+ Ckij^f + ^(f ijg ij ) - 12vrf o - 24^V + ^(%<W 
Using the constraint (|4.2j) we get 

^M S ) = (^HM^H^Too+Svr^^ (4.15) 

Thus using (|4.8|) and the definitions of e and S, we obtain 

d_ 

dr 

The last step is due to the fact that considering expanding spacetimes implies 
(kijg l i) < 0. In |2j it is shown that there exists a lower bound of e, say e*, 
which only depends on the initial condition of the spacetimes. Therefore we 
have 

S(r) = 0(e-^) 
and as a consequence from the definition of S we have 

(k l3 9 ij ) 2 - 12n(ij 2 + 2Vo) = 0(e~ € ' T ) (4.16) 

By (gHH) the rest of the claims follows ; a^a^ = R = f 00 = 0(e~ t * T ). Here 
while we use the notation O(-), we lose a piece of information from above 
that S is non-negative. So we want to point out that the errors in ()4.9(l - 
(|4.12|) are non-negative. □ 



-S = -eS + 2(k ij g ij )[(a ij a i i) + 4m{T ijg ^) + 47rf 00 ] < -eS 
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4.1. Asymptotic behaviours of (k^g 13 ) and t/j. The estimate (|4.1(ijl is 

unsatisfactory in the sense that we do not have sufficient information to say 
individual asymptotic behaviours of (kijg 1 - 1 ) and ■0. In this subsection, we 
obtain asymptotic behaviours of these quantities, separately. 

Proposition 4. 

= A 7 + 0(e-v T ) (4.17) 
k ijg v = -3k 7 + 0(e- VT ) (4.18) 

where 7 := v / 2To/v6 — A 2 and 77 := min{i«;7(6 — A 2 ), e*/2} 

Note that in the case < A < a/2/3, we have 77 = e*/2. 1 
Proof : With the evolution equation (|4.6|) consider 

where 

- A 7 ) := (fcy<7« + 3k 7 )(^ - A 7 ) + y$ ~ X ^ 

+ (A 2 - 3)k 7 (0 - A 7 ) + \j(kj9 ij + 3k 7 ) 

(|4.9j) implies 

{h^f = k 2 - A7) 2 + 3A 7 (^ - A 7 ) + 9 7 2 ] + 0(e" eV ) 

When = A7, we get 

+ 3k 7 = 3k 7 - a/9k 2 7 2 + 0(e~ e * r ) = 0(e~ e * T ) 

Then -F^(O) = 0(e~ e * T ) and after a lengthy elementary computation one 
can see that ^(0) = —7^7(6 — A 2 ) + 0(e~ e * T ). So when is close to A7, 

A 7 ) = 0(e~^) + [-i« 7 (6 - A 2 ) + O(e" eV )](0 - A 7 ) 
+ C§ - A 7 ) 2 
Define Y$(t) := e^^ 6 "^) 7 "(^ - A7). Then 



e*, the lower bound of e depends on not only the constants A and Vo but also initial 
data. When A £ (0,^/2/3), the trivial lower bound of e, which may not be sharp, is 
iK V / 6V r (2 - 3A 2 ) (see [HOI for details). In this case ±^7(6 - A 2 ) > e*/2 is true. 
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This yields 
This implies that 

Consequently, 

5^ = C + 0{e~^- x2)T ) + e^ K7(6 - A2)T 0(e- e * T/2 ) 
Therefore we have 

^ _ A 7 = Ce-5 K 7(6-A 2 )r + ^^-^(e-A^r) + 0{e~ e * T ' 2 ) 

and this gives the proof of (|4.17|) . For (|4.18f) , with (|4,8() and (|4.15f) consider 
^(hj9 ij + 3k 7 ) = FjiM* + 3k 7 ) 

where 

■= \{k l9 13 + 3 K7 ) 2 + - A 7 ) 2 + ^(kjgv + 3k 7 )(^ - A 7 ) 

+ ^X 7 k 2 (^ - A 7 ) + i(A 2 - 4)k 7 (%^' + 3k 7 ) 
+ (a ij a ij )+4nT 00 + 4Tr(f ij g i i) 

implies 

V" 2 = ^2^9 ij + 3^7) 2 " ^(kijg ij + 3k 7 ) + A 2 7 2 + 0(e^* T ) 

When kijg 1 ^ = — 3k 7 , one can see that F^(0) = 0(e~ e * T ) and also F^(0) = 
-|k 7 (6 - A 2 ) + 0(e" e * T ). Therefore 

^(%<? ij + 3*7) = 0(e~^ T ) + [-~kA(6 - A 2 ) + 0( e - eV )](%^' + 3k 7 ) 
+ + 3* 7 ) 2 

where kijg lJ +3k 7 is small and by the same argument as above (|4.18f) follows. 
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4.2. Relation between r and t and asymptotics in terms of t. So 

far in the present section, we have obtain asymptotics of quantities, a^a 13 , 
R, Too, kij and gij, in terms of the time coordinate r after rescaled by a 
certain factor of the scalar field. In order to study further asymptotics, it is 
necessary to recover these quantities in terms of the time coordinate t. For 
this reason, in this subsection the relation between the two time coordinates 
and the rescaling factor g - ^^/ 2 in terms or t will be analysed. 



Proposition 5. 



-A 2 K7t/2 



2 e AftCi/2 

A 2 k 7 



0(t- 2 lnt), if\ 2 K~f/2 = r) 
0(t-^+V), if\ 2 n~f/2>ri 
{ 0(t- 2 ), if\ 2 K 1 /2<i 1 



where Q := 27//A 2 k 7 and C% is a constant. 

PROOF : From (|4.17|) we have 

<f>(r) = A 7 r + C x + 0{e-^ T ) 

So by (|4~TT) we get 



t(r) = t(T ) + e 
where t(jo) = io- Then 



AkCi/2 



exp[A 2 K 7 s/2 + 0(e- r ' s )]ds 



-A 2 K7r/2 



*(t) 



'X 2 K'Jt/2, 



AkCi/2 -A 2 k7t/2 



-A 2 K7t/2 



t(r ) + e 



AkCi/2 



/ ex P 

[A 2 K7s/2 + C(e" 
1 + 0( e -" s )] ds 



'To 

= A 2 K7 ( S -r)/2 



TO 



2e A K C*i/2 



+ < 



A 2 k 7 

In all cases, we have 
Consequently 



' 0(Te~ x2 ^ T / 2 ), ifA 2 K 7 /2 = r/ 
0(e~^ T ), HX 2 K-//2>7] 
k 0(e- A2 ^ r / 2 ), if A 2 k 7 /2 < 77 

e -A 2 K7 r/2^ T ) < c 



-A 2 K7r/2 



or 1 ] 



(4.19) 



Thus the proposition follows. 



□ 
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Proposition 6. 

0(f- 2 lni), if\ 2 Kj/2 = ri 

0(f-^ +1 )), if\ 2 Kj/2>r] 

0(t- 2 ), if\ 2 K 7 /2< V 

Proof : By KWi . 

e -\K<t>/2 _ e -A 2 K 7 r/2 e -AKCi/2^ + Qf e -rrr\\ 

Combining this with Proposition yields the conclusion of the proposition. 

□ 

Proposition 7. 

aij a ij = 0(r (5+2) ) (4.20) 

R = 0{t-^ +2) ) (4.21) 

Too = 0{t-^+V) (4.22) 

where £ := 2e* /\ 2 k^. 

Proof : By Proposition El 

e -"* T = o{t- 2e *' x2 ^ T ) 

Also Proposition |S] implies 

Combining these with (|4.1()j) - (j4.12j) in Proposition 01 concludes the propo- 
sition. □ 



-\ncf>/2 



4.3. Asymptotic behaviours of kijg^ , tp and ^ in terms of t. In this 
part, we will observe asymptotic behaviours of kijg^ , ip and 4> in terms of t 
using the relation between two time coordinates r and t we have obtained 
in the previous subsection. 
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Proposition 8. 



u. .Jj 



Xk 



6 l 
Y2 t H< 



0(t~ 2 lnt), 
0(t- 2 lnt), 

o(Hf +1 )), 

l or 2 ), 



if X 2 kj/2 = rj 
if X 2 kj/2 > 7] 
if X 2 kj/2 < 7] 

if X 2 K~f/2 = r] 
if X 2 kj/2 > rj 
if X 2 kj/2 < rj 



(4.23) 



(4.24) 



Proof : Note that from Proposition [3 

So combining ()4. 17|) in Proposition|I]and Proposition|Bl (|4.23|) follows. With 

□ 



1)4.18(1 in Proposition 0] the same argument applies to prove 1)4.24(1 
Proposition 9. 



Xk 



lnt + C 2 + < 



0(t-Hnt), 

o(t- 1 ), 



if X 2 K'y/2 = rj 
if X 2 kj/2 > rj 
if X 2 kj/2 < rj 



where C2 is a constant. 



Proof : The proposition follows directly from 1)4.23)) in the previous 
proposition. □ 

4.4. Asymptotic behaviours of g- L j and g %0 . Asymptotics of gij and g 13 
will be analysed first in the time coordinate r. By means of the relation be- 
tween the two time coordinates r and t in Subsection l4.21 these asymptotics 
will be recovered in terms of t. 

It has been observed in Proposition |3] that a^a 13 = 0(e~ e r ). Using the 
following two lemmas we will identify &ij{r) which play a role to analyse 
9ij{r). 

Lemma 1. Let ||<?(t)||, ||fc(T)|| and \\ct(t)\\ denote the norms of the matrices 
with entries gij(r), kij(r) and aij(r), respectively. Then 

\\a(r)\\<Ce-^/ 2 \\g(r)\\ 

Proof : The lemma follows by the fact that 

[|*(r)[U T) <(^) 1 / 2 <Ce- e,T / a 
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The last inequality is due to (|4.1()j) . □ 
Lemma 2. 

|e- 2K7T ^(r)| < C (4.25) 
\e 2K ~ fT g ij (T)\ < C (4.26) 



Proof : Let gij(r) := e~ 2K7T ^(r). Then using (|4"%|) and (|1~Tg)) . we get 

— ^ = -2k 7 ^ - -{hm9 lm )9ij ~ 2e" 2re7T ffi j 

= 0(e"" T )^-2e- 2K ^ (4.27) 

Now let us use the norms again. Let ||<7(t)|| be a norm of the matrix with 
entries gij(r). Then we get 

||fif(r)[| < \\9(ro)\\+ [ T (Ce-^\\g(s)\\+Ce~ 2 ^ s \\a(s)\\)ds (4.28) 

J To 

< \\9(ro)\\+ r(Ce-^||5( S )||+C e - 2 ^||a( S )|| § ||g( S )||)(i S 



TO 



Note that 

P(t)\\~9 < e 2 ^ T (aija ij ) 1/2 < Ce 2 ^ T e~ € * T/2 
So combining this with (|4.28j) yields 



||ff(r)||<||5(ro)||+ / Ce^ s \\g(s)\\ds 
J ■ 

By Gronwall's inequality, this becomes 



TO 



\g(r)\\ < ||g(r )||exp 



Ce-i" ds 

TO 



< c 



Therefore |ffij(r)| is bounded by a constant for all t > tq. Also (|4.26j) follows 
by the same argument. □ 

Proposition 10. 

e-^aa = 0{e- e ' T/2 ) 

Proof : The proposition follows by Lemmas ^ and El □ 
Proposition 11. 

g ij (T) = e 2K ^(g ij + 0(e-^)) (4.29) 

& «(r) = e- 2 ^ T (g ij + 0(e-^)) (4.30) 
Here Qij and Q l i are independent of r. 
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PROOF : Let us consider (|4.27j) again ; 

^•(r) = 0(e-^)g t ,(r) - 2 e - 2K ^(r) 

where gij(r) = e~ 2K ^ T gij(r). Then Lemma El and Proposition 1101 imply 

^Qijir) < Ce-^gij + Ce~^ r / 2 < Ce~^ T 

Since 4f(jij is decaying exponentially, there exists a limit, say Gij , of g^ as 
r goes to infinity. Then this gives 

~9i j {r)=G ij + 0(e-^) 

i.e. 

Here the lower order term of gij(r) is of an exponential form so that it is 
combined with the leading order term, which makes it possible to compute 
explicitly. So <? u (t) is 

0«( r ) = e~ 2 ^ T {g ij + 0(e-" T )) 



□ 



Proposition 12. 



9ij 



(t) = t 4 ^ 



X 2 KJ \4A : 



V 



2 6 AkCi/2 



Sin + < 



Oft- 1 hit) 

I or 1 ), 



if X 2 K r y/2 = r\ 
if \ 2 kj/2 > rj 
if X 2 kj/2 < 7] 



(4.31) 



/2p AkCi/2 4 / A 2 

(W) S *' + < 



' Oft -1 bit), if\ 2 Kj/2 = r) 



I or 1 ), 



if X 2 K-f/2 > rj 
if X 2 kj/2 < rj 



(4.32) 



Proof : Recall that 

e -vr = ^~C) 
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Proposition 03 implies 



4/A 2 



+ < 



V 



V 



2gAftCi/2 

A 2 k 7 

A 2 K 7 \ i A J 
2 e A«;Ci/2 



' ©(rMnt), ifA 2 K 7 /2 = r ? 

0(rC), ifA 2 K 7 /2>r/ 
_ ©(r 1 ), if A 2 k 7 /2 < 77 

' Oft" 1 In i), ifA 2 K 7 /2 = T/ 
+ < O(t^), ifA 2 K 7 /2>?7 
k Oit' 1 ), if A 2 k 7 /2 < 77 



So by means of this and Proposition II H 1)4, 31ft follows. The same argument 
for gV yields (JQ2l . □ 

4.5. Asymptotics of the generalized Kasner exponents and the de- 
celeration parameter. Let Aj be the eigenvalues of with respect to 
gij(t), i.e., the solutions of det(kj — X6j) = 0. Define the generalized Kasner 
exponents by 

_ Aj _ A, 

ft : " " (W m ) 

The name comes from the fact that in the special case of the Kasner solutions 
these are the Kasner exponents. Note that while the Kasner exponents are 
constants, the generalized Kasner exponents are in general functions of t. 
The generalized Kasner exponents always satisfy the first of the two Kasner 
relations, but in general do not satisfy the second, where these two Kasner 
relations are 



52 Pi 

i 



1. 



(4.33) 
(4.34) 



The following proposition exhibits that the spacetime isotropizes at late 
times. 



Proposition 13. 



where £ = 2e*/A 2 K 7 . 



Pi(t) = - + 0(i 



Proof : First note that by 1)4 . 8 j) Aj are also the solutions of 



det {a)-[\--{k lm g lm )]8f ) =0 
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So the eigenvalues of o"y(r) with respect to gij(r) are 

A, := Xi - \{h m g lm ) 

Also note that Y.i&i? = oi m o lm '. Then (ETTUl) implies 

Therefore using this and (|4.18f) we obtain 

„.-- = ^ = 0(e- e * T/2 ) 

Thus Proposition [5] completes the proof. □ 

There is another quantity to be considered regarding expanding cosmo- 
logical models, which is the deceleration parameter, say q. This deceleration 
parameter is related to the mean curvature, as follows 

d_ 

In accelerated expanding models, the deceleration parameter is negative. 



Proposition 14. 

' Oft- 1 hit), if \ 2 Kj/2 = r) 

0(t-<), if\ 2 K7 /2> V 

k Oit- 1 ), if\ 2 K 1 /2<i 1 



X 2 



Proof : The proof is a straight forward computation from (|4.24l) . □ 

4.6. Asymptotics of momenta. We will analyse the behaviour of the 
momenta of the distribution function / along the characteristics where / 
is a constant. From the Vlasov equation (|2.11j) we define the characteristic 
curve V i (t) by 

— = 2k)Vi - (1 + g rs V r V s r 1/2 lLv m V n (4.35) 

for each V l (to) = Vq given to. The characteristics Vi, rather than V 1 , has a 
simpler form, which makes analysing the behaviour of the momenta easier. 
So here Vi(t) satisfies 

^ = -(1 + g rs V r VT 1/2 lLnV P V q9 P m g^ 9l , (4.36) 
for each Vi(to) = viq given to- Also observe that Vi(r) satisfies 

^ = - e ^/ 2 (l + g rs v r v s )- 1/2 YmnypV q 9 pm g gn g ll (4.37) 
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For the rest of the paper, the capital V % and Vi indicate that v % and vi are 
parameterized by the coordinate time t or r. 



Theorem 3. Vi(t) from \4-Sb\) converges to a constant along the character- 
istics as t goes to infinity. That is 

Vi{t) - C ' i+ \0 ( t- % ,/ C >l, 
where u := min{£,4/A 2 — !}■ Furthermore, 



V*(t) = t 



-4/A 2 



/ f Oit-Hnt), if\ 2 Kj/2 = r] 
C 4 + I 0(t-<), if\ 2 ^/2 > r, 

\ [o(t~ l ), if\ 2 K 7 /2<n 



Before the proof of this theorem, some lemmas are required. 
Lemma 3. 

{g^ViV^r) = e~ 2 ^(V + C(e-^)) 
where V is a constant. 

Proof : First note that by Propositions I1UI and 111! we have 

Let a 13 := e ( 2K7+e * l 2 ) T a %3 . Then a 13 is bounded by a constant for all r > r . 
Since Q %3 in Proposition ^2 is positive definite and time independent, there 
exists a constant C, independent of time, such that 

ff U'.ij < (V' j \)Vj 

Then by means of this and we obtain 

d 

~aW' 

2 

= 3< 

< ( - 2k 7 + Ce - " 1 ")^'^ + Ce-^^'^G^ViVj 

< ( - 2k 7 + Ce"^) ^Vi^- (4.38) 

Here to get the first equal sign, ()4.37|) is used. Yet the terms involved with 
(|4.37|) vanish due to the antisymmetric property of r j\ nn combining with g %3 . 
Now consider V T := e 2K ~< T g ij ViVj . Then one can see from (J4~38|) that 

dr 



-{'i'Y.Yn = 2k r <Y,Y, 



;(k m g lm )9 ij V l V j + 2a i iV i V 3 
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So there exists a limit of V T , say V as r goes to infinity. Then we have 

V T = (V + 0{e-^)) 
This completes the proof. □ 
Lemma 4. 



/ 2p \kC 1 /2\ 4/A " 

(9"ViVj){t) = V \— ? —\ t 



4/A 2 



C(t- 2 lnt), if\ 2 Kj/2 = r] 
+ < i/A 2 K 7 /2 > 77 

k 0(t- 2 ), z/A 2 K 7 /2<r? 



Proof : Recall that 



-77T 



0(f 



Then by means of Proposition and Lemma El we have 

4/A 2 



G^V^-Xt) = v 



2e A«;Ci/2 



A 2 



+ o(t 



-(c+A) 



+ < 



0(r 2 lnt), ifA 2 K7/2 = r? 
0(rK +1 >), ifA 2 K 7 /2>r/ 
[ 0(r 2 ), if A 2 k 7 /2 < 7? 



□ 



Note that C + 4/A 2 > ( + 2. So the lemma follows. 
Lemma 5. 

\vm <c 5 + o(t-<) 

where C§ is a constant and for all i. 

Proof : Since Q tJ is positive definite, there exists a constant C such that 

W 2 (T)<C^Vi(r)^(r) 

Note that ( 2f ~ : C/J\)\- j [ s the leading order term in (g ij ViVj)(T) in (jO0)l . 
So using Lemma El we conclude that 

g^V^- = V + 0(e-* r ) 

Combining this with the fact that 

e-^ = 0(r<) 

we complete the proof. □ 
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Proof of Theorem El Note that in Bianchi type I, since all structure 
constants are zero, also the Ricci rotation coefficients 7m m are zero. So from 
(|4.36|) it is clear that Vi(t) = for all t. 

More generally Lemma El says that all V* are bounded by a constant when 
t goes to infinity. However this allows oscillating behaviours. So to rule out 
these cases, it is necessary to analyse -jf. Combining Proposition 1121 and 
(|4.3fi|) . we have 



dVi 



dt 



< ct~ A ' x2 \v p \\v q \ 



Here the right hand side is a summation for some p and q. Then by Lemma 
[S] this implies 

dV t 



dt 



(4.39) 



This leads to the conclusion that when t goes to the infinity, ^4- goes to 
zero, and so Vi goes to a constant. Combining Lemma |S] and Q4.39JI we have 



Vi(t) =c 3 + o(t-<) + o(t 



-(4/A 2 -^ 



Since 4/A 2 



1 > 1 



Vi(t) =c 3 + 



o(t~ 



if C < 1, 
if C > 1, 



where lo =: min{4/A 2 — 1, C}- Now combining this with (|4.32|) completes the 
proof. □ 



4.7. Geodesic completeness. In this part, we will prove the completeness 
of future directed causal geodesies, which has been postponed in Subsection 

Proof of Theorem 21 The geodesic equations for a metric of the form 
(|2.1j) imply that along the geodesies the variables t, v % and v° satisfy the 
following system of differential equations : 

dt 
ds 
dv° 
ds 
dt/ 
ds 



(4.40) 



kijV l v 3 



^- = 2kyv°-Y mn v m v n 
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where s is an affine parameter. For a particle with rest mass m moving 
forward in time, v can be expressed by the remaining variables, 

v° = (m 2 + <^V) 1/2 (4.41) 

The geodesic completeness is decided by looking at the relation between t 
and the affine parameter s, along any future directed causal geodesic. This 
relation is clear from (|4,4()|) and (|4.41jl . I.e., it is given by 

| = ( m 2 + ^ V) l/2 

To control this, it is necessary to control gijV l v 3 as a function of the coor- 
dinate time t. Consider first the case of a timelike geodesic. I.e., m > 0. 
Then V*(t) satisfy 

In this case, the arguments presented in Subsection 14 . 61 is valid, in particular 
those in Lemmas 01 and 0] when m = 1. Therefore by LemmaEl (gijy l V :, )(t) 
is bounded above by Ci~ 4 / A2 , and so by C, for all t > to. Hence this gives 
that 

f t = (m 2 + ( giJ V^)(t))- l/2 >C 

Therefore when s is recovered by integrating this, the integral of the right 
hand side diverges as t goes to infinity. 

Now consider a null geodesic, i.e., m = 0. Then in this case V l {t) satisfy 
dV i 



dt 



2k)Vi - (gijV^r^^V^ 



Also Lemma 0] is valid. Therefore (gijV l V J )(t) is bounded by a constant 
and this gives 

% = {m{v l v 3 ){t))- 1/2 >c 

Therefore as t goes to infinity so does s. □ 



4.8. Asymptotics of the energy-momentum tensor. In this subsection 
it will be analysed the asymptotic behaviour of the energy-momentum tensor 
in an orthonormal frame on the hyper surf aces. 
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Proposition 15. Let{ii} be an orthonormal frame. The energy-momentum 
tensor is described by 

p(t) = J f(t,v)(l + \v\ 2 )V 2 dv (4.42) 
Jiit) = j f(t,v)vidv (4.43) 
Sij(t) = J + \v\ 2 )' 1/2 dv (4.44) 

where p := Tqq is the energy density, Ji := Tq{ the components of the current 
density and S{j := Tij are the spatial components of the energy-momentum 
tensor. Here the hats indicate that objects are written in the orthonormal 
frame. Furthermore v := (vi,V2,va) and dv = dvxd^dv^. 

Then pit), Jiit) and Sij(t) tend to zero as t goes to infinity. More pre- 
cisely, 

P (t) = o(t-^ 2 ) 
Ji(t) = o(t- 8 / x2 ) 
s tJ (t) = o(t- 10 / x2 ) 



Furthermore 



p(t) 



0{t' z ' x ) (4.45) 



^# = 0(r 4 / A2 ) (4.46) 
pit) y ' K J 

Proof : Note that f(to,v) has compact support on v. Also observe that 
Theorem |31 implies that V{(t) is uniformly bounded. Combining these two 
facts implies that there exists a constant C such that 

f(t,v)=0, if H>C (4.47) 

for all t. By (jOTj) we have 

f(t,v) = 0, if \vi\>Cr 2/x2 (4.48) 

So using ()4.47j) and (|4.48|) we get 

p(t)= I fit^il + ^^dv 

J\v t \<Ct- 2 /* 2 

Note that since f(t,v) is a constant along the characteristics, 

< ll/oll := sup{|/(i 0j i5)| : for all v} (4.49) 
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So we obtain 



p(t)<C f(t,v)dv 

>\vi\<Ct- 2 /* 2 



< C||/ p- 6/A2 
Also by (|Q7|) - (pQ9|) we have 



Ji(t) = I f{t,v)vidv 

< Ct~ 2 / x2 [ f{t,v)dv 

J\v,\<Ct- 2 / x2 

<c\\f \\t- 8 / x2 



Similarly 



Sij(t) = / f(t,v)viVj(l + \v\ 2 r 1/2 dv 

J\vi\<Ct- 2 />> 2 

< Ct~ 4 / x2 [ f{t,v)dv 

J^KCt- 2 ^ 2 

<c\\f Q \\t~ w ' x2 

Now let us estimate the ratios Ji/p and Sij/p. By means of (|4.47j) and 
(g3HJ), we get 

Jj{t) _ J f(t,v)vj dv 
Pit) ~ J f(t,v)(l + \v\ 2 y/ 2 dv 

< rt -2/x 2 ff(t,v)dv 

j f(t,v)(i + \v\ 2 y/ 2 dv 

< ct- 2 ' x2 



Similarly 



P {t) j f{t,v)(i + \v\ 2 y/ 2 dv 
< ct- 4 ^ 2 



□ 



In this proposition since all components of the energy momentum tensor 
in an orthonormal frame go to zero as t goes to infinity, it can be concluded 
that in a certain sense solutions of Einstein- Vlasov system coupled to a non- 
linear scalar field with a exponential potential are approximated by vacuum 
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Einstein solutions. In a more detailed level (j4.45|) and (|4.46|) resemble the 
non-tilted dust-like solutions in which Ji{t) and Sij(t) are identically zero. 
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